fig = figure;
errorbar(trendp,M50(:,3),M50(:,3)-M50(:,1),M50(:,5)-M50(:,3),'LineStyle','none','LineWidth',1,'Color',rgb('LightGray'));
hold on;
  errorbar(trendp,M50(:,3),M50(:,3)-M50(:,2),M50(:,4)-M50(:,3),'LineStyle','none','LineWidth',1,'Color',rgb('Black'));
hold off;
xlim([0 114]);
ylim([-0.02 0.09]);
orient landscape;
print(fig, [figdir, 'Figure_5_a_50yr' mod_name], '-dpdf','-bestfit');
close(gcf);

fig = figure;
errorbar(trendp,M100(:,3),M100(:,3)-M100(:,1),M100(:,5)-M100(:,3),'LineStyle','none','LineWidth',1,'Color',rgb('LightGray'));
hold on;
  errorbar(trendp,M100(:,3),M100(:,3)-M100(:,2),M100(:,4)-M100(:,3),'LineStyle','none','LineWidth',1,'Color',rgb('Black'));
hold off;
xlim([0 114]);
ylim([-0.02 0.09]);
orient landscape;
print(fig, [figdir, 'Figure_5_b_100yr' mod_name], '-dpdf','-bestfit');
close(gcf);